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Abstract 

The electroweak form factors of the nucleon as obtained within a three flavor 
pseudoscalar vector meson soliton model are employed to predict the ratio of the 
proton and neutron yields from 12 C, which are induced by quasi-elastic neutrino 
reactions. These predictions are found to vary only moderately in the parameter 
space allowed by the model. The antineutrino flux of the up-coming experiment 
determining this ratio was previously overestimated. The corresponding correction 
is shown to have only a small effect on the predicted ratio. However, it is found that 
the experimental result for the ratio crucially depends on an accurate measurement 
of the energy of the knocked out nucleon. 







1. Introduction 



By means of a continuum random phase approximation (RPA) JT|, || it has been shown 
0, fl] that the ratio R y of the proton-to-neutron neutrino-induced quasi-elastic yields from 
an isoscalar nucleus provides a relation between the nucleon matrix elements of the strange 
quark vector and axial vector currents (A^|s7 M s|A^) and (-/V|s7 M 7 5 s|./V), respectively. This 
is especially interesting because it allows one to determine the singlet matrix element 
of the axial vector current (the so-called proton spin puzzle ||) without assuming flavor 
SU(3) symmetry This symmetry is commonly employed to relate matrix elements of 
currents which conserve the flavor charges to those which change this quantum number. 
The latter are measured in semi-leptonic hyperon decays. 

The purpose of the present paper is three-fold. Firstly, we will update the earlier 
calculations || |4|] by substituting more accurate energy spectra for the LAMPF neutrino 
and antineutrino beam which recently became available 0. This spectrum enters the 
computation of the quasi elastic yields by folding it into the cross-section for the (anti) 
neutrino - nucleus reaction 0]. In the second part we will adopt the SU(3) Skyrme model 
with vector mesons to compute the relevant momentum dependent form factors of the 
nucleon. This makes possible a prediction for R y because these form factors serve as the 
essential ingredients for the nuclear model calculation developed in refs. |3], |J. Although 
the incorporation of vector meson fields into the SU(3) Skyrme model represents a severe 
complication || these fields are unavoidable for a realistic description of static baryon 
properties. In particular these fields account for short range effects, which are particularly 
important to gain a realistic, non-zero value for ||. The simple Skyrme model of 
pseudoscalars predicts gA{<f) = |]TIJ. Such short range components are also necessary to 
gain a non-zero result for the neutron-proton mass difference [|llj . On the basis that this 
model reproduces the experimental two-flavor form factors only within an error of 20- 
30% (which already represents an improvement compared with the pure SU(3) Skyrme 
model) one might want to argue that other approaches may be more appropriate to 
compute strange form factors. However, this model is unique in the sense that it treats 
vector and axial vector currents simultaneously and contains a consistent and systematic 
treatment of symmetry breaking effects. The resulting predictions for R y turn out to be 
quite insensitive to the model parameters. This indicates that the strange form factors 
may independently be considered from those in the two flavor subspace. A recent study 
[p~3]| has shown that the ratio R y is not very sensitive to the nuclear model but significantly 
depends on the prediction for (N \ s\N). In this respect a precise measurement of R y 
by the LSND-Collaboration at LAMPF [13 will be discriminating on the nucleon model. 
In section 3 we will briefly review the computation of the relevant form factors in this 
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model after having outlined (in section 2) their relevance for the (anti)neutrino scattering 
off nuclei in the context of the continuum RPA. 

Finally we will point out a possible source of uncertainty in the experiment associated 
with the sensitive dependence of R y on the energy thresholds for the integrated proton 
and neutron yields. The corresponding results together with a comprehensive discussion 
will be presented in section 4. 

2. The u — if nuclei reaction 

At low energies (compared to the electro-weak scale) the interaction of neutrinos with 
matter is described by the current-current interaction introducing the Fermi coupling 
constant Gf 

4 ff = -|(^ ( - ) +jiv (o) +h.c.). (i) 

Here denotes the leptonic part of the electro-weak current^] while J^w represents 
the hadronic counterpart. Of course, these currents break parity in a maximal way, i.e. 
jjfl) = v (a) _ a (a) wr th v anc [ a referring to the vector and axial-vector currents. The 
hadronic (axial)vector currents (V and A) are defined similarly. 

To describe neutrino scattering off nuclei in a perturbative treatment of £ e fr we will 
employ the continuum RPA approach (as a mean field model) of ref.JT], [2| in order to 
compute matrix elements of when sandwiched between nuclei states^. In this nuclear 
model the interaction between the constituents of the nucleus is described by combining 
the usual RPA treatment with a correct description of the particle states in the continuum, 
i.e. the excited many-body states are coherent superpositions of one-particle-one-hole (lp- 
lh) excitations obeying the proper Coulomb boundary conditions for scattering states. 
For details on this approach we refer to the literature [I], ||. Here we restrict ourselves to 
a brief description of its basic features and ingredients. The continuum RPA provides a 
good description of the nuclear ground state while the excited states are generic continuum 
states possessing a lp-lh structure. Final state interactions are accounted for by a realistic 
(finite range) residual interaction derived from the Bonn meson exchange potential |T5| , 
It should also be mentioned that this treatment assumes that the reactions proceed 
predominantly via the knocked out nucleon. As the main result of the continuum RPA 
approach the scattering cross section of the (anti)neutrino nuclei reaction is expressed in 
terms of the nucleon matrix elements 

(N\J^ \N) , a = 0,..,8, (2) 

1 This shorthand notation is assumed to include the dependence on the Weinberg angle, 6\v- 
2 The expression for the scattering cross section is given in eq (11) of ref. [Q. 
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where \N) represents a nucleoli isospinor. Since the nucleon is an extended object form 
factors for this matrix elements are mandatory. The Lorentz covariant expressions for the 
(axial)vector current matrix elements^] of the flavor conserving components (a = 0, 3, 8) 
read 



(N'\VZ\N) = u N (p') \f? (V) 7m + Ft (V) ^ 



u N (p) 



^2M_ 

(N'\A«\N) = u N (p') [G a A (q 2 ) 7/ , 75 ] u N (p) (3) 

where u^(p) denotes the proton Dirac spinor. The neutron form factors are obtained from 
isospin covariance. The linear combinations relevant for the processes under consideration 
are @, § 

F t (q 2 ) = ± Q - sin 2 6 w ) Ff{q 2 ) - sin 2 e w i?=V) - ^(q 2 ) (i = l,2), 

G(q 2 ) = T\G\{q 2 )+ l -G\{q% (4) 

where 1 = and s denote the non-strange and strange combinations of a = 0,8, respec- 
tively. The upper (lower) sign refers to the proton (neutron). The vector form factors Fi 
are more conveniently expressed in terms of electric and magnetic combinations 

G a E {q 2 ) = F«{q 2 ) - j^F^q 2 ) and G a M (q 2 ) = F?{q 2 ) + F«(q 2 ). (5) 

The charge matrix is given by Q = (A 3 + A 8 /v / 3)/2 hence one may adopt the physically 
more transparent basis consisting of the electro-magnetic form factors of the nucleon, 
G^ n M (q 2 ). Due to isospin covariance G\(q 2 ) is identical to the axial charge GA(q 2 ) of the 
nucleon extracted from neutron /3-decay. The computation of these form factors in the 
Skyrme model with vector mesons will be reviewed in the next section. 

3. Nucleon form factors in the vector meson Skyrme model 

Our starting point is a three-flavor chirally invariant theory for pseudoscalar and vector 
mesons. The model Lagrangian also contains abnormal parity terms to accommodate 
processes like u —>■ 3tt. A minimal set of symmetry breaking terms, which transform as 
(3 x 3 + h.c.) under the chiral group SUl(3) x SUr(3), is included [TIJ to account for 
different masses and decay constants. This effective theory contains topologically non- 
trivial static solutions, which are constructed by imposing ansdtze in the isospin subgroup 

fi \ G(r) 

£„(r) = exp ( -r ■ rF(r)J , u (r) = u{r) and p i>a (r) = —^—^ij a rj , (6) 



3 We omit the induced form factor of the axial current. 
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while all other field components vanish. Here ^ = exp (iir ■ T/2f n ) refers to the non- 
linear realization of the pion fields. The radial functions are determined by extremizing 
the static energy functional subject to boundary conditions appropriate to the topolog- 
ical sector of winding number one. Motivated by the large Nc studies of QCD these 
solitons are identified as the baryons with unit baryon number [O. Unfortunately the 



field configuration (|B|) carries neither good spin nor flavor quantum numbers, hence an 
appropriate projection has to be performed. Also strange flavors are not yet excited. 
These two features are accounted for by introducing time dependent collective coordi- 
nates A(t) for the (approximate) zero modes associated with the SU(3) vector rotations: 
£(r,t) = A(t)£ n (r)A^(t) and similarly for the vector meson nonet. Furthermore field 
components, which vanish classically like e.g. the kaons, are induced by the collective 
rotations. These induced fields turn out to be proportional to the angular velocities 

i 8 

A^A=-Y, ^ a (7) 

^ a =l 

and exist for both, rotations into strange (a = 4, .., 7) and non-strange (a = 1, 2, 3) direc- 
tions. In order to compute the induced fields ansdtze are chosen, which are consistent with 
the parity as well as Lorentz and isospin structure of the considered mesons. Altogether 
eight real radial functions are introduced, which solve inhomogeneous linear differential 
equations derived from a variational principle to the moments of inertia M. The classical 
fields (P) act as source terms, which stem from the abnormal parity parts of the effective 
meson theory. The collective Lagrangian L = L(A, {Q a }) is extracted and the collective 
coordinates are quantized canonically 

dL(A, {&}) 

Ra = M~« (8) 

providing a linear relation between the angular velocities Q a and the right generators of 
SU(3), R a . The fact that these objects are operators is reflected by the commutation 
relations^] [R a ,Rb] = —f a bcR c , with f abc being the structure constants of SU(3). The 
resulting collective Hamiltonian may be diagonalized exactly [19| yielding the spectrum 
of the low-lying | + and | + baryons. In the absence of flavor symmetry breaking the 
eigenf unctions *&(A) of this Hamiltonian are SU(3) D-functions associated with a certain 
representation (e.g. the 8 for the nucleon). Once flavor symmetry breaking is included 
these D-functions become distorted reflecting the admixture of higher dimensional repre- 
sentations as for example the 10 or 27. 

Extending the action to account for electro-weak interactions allows one to derive 
covariant expressions for the (axial-)vector currents from the terms linear in the cor- 

4 The right generators of SU (3) are identified in the body fixed system of the rotating soliton. 
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Table I: Predicted form factors at zero momentum transfer. Three sets of parameters 
in the model Lagrangian are used, (x): best fit to the baryon mass differences, (*): 
reproducing the experimental value for Ga(0), (+): best fit to Ga(0), G p m (0) and G^j(O). 
Furthermore the results associated with 577(3) symmetric wave-functions are presented 
for the parameter set (x). 





G A (0) 


G P M (0) 




G A (0) 


G S M (0) 


X 


0.94 


2.35 


-1.86 


-0.030 


-0.055 


* 


1.25 


3.23 


-2.86 


-0.015 


-0.035 


+ 


1.02 


2.57 


-2.11 


-0.028 


-0.050 


SU(3) sym. 


0.88 


2.48 


-1.54 


-0.058 


-0.559 


Exp. 


1.25 


2.79 


-1.91 







responding gauge fields. Substituting the above described ansatze and eliminating the 
angular velocities in favor of the generators ([8]) leaves the currents as linear combinations 
of radial functions and operators in the space of the collective coordinates A. The former 
are given in terms of the classical and induced profile functions of the meson fields. The 
matrix elements of the currents, which eventually provide the relevant form factors, are 
computed in two steps. Firstly, the radial functions are Fourier-transformed yielding the 
momentum dependence of the form factors in the Breit frame |[20|| . Secondly, an SU(3) 
"Euler-angle" representation for A(t) is employed [|19[] to parametrize the SU(3) operators 
as well as the exact eigenf unctions ^{A). This makes possible the evaluation of the spin 
and flavor parts in the matrix elements of the currents. 

In ref. |§ the parameters of the model Lagrangian, which could not be determined 
from the meson sector, were adjusted to provide a best fit to the baryon mass differences. 
This has lead to reasonable, though not perfect, agreement for the electro-magnetic and 
axial form factors. The strange form factors, which are obtained by considering different 
flavor components of the currents, were predicted. In table | we display the results for 
zero momentum transfer together with results corresponding to parameter sets which 
improve on the axial charge of the nucleon and/or the magnetic magnetic moment of 
the nucleon. Also displayed are the predicted strange form factors at zero momentum 
transfer. Generally we find that all parametrizations predict the strangeness form factors 
to be negative and small in magnitude. In fig. [I] we display the predicted momentum 
dependence of the strange vector form factor Gs{q 2 )- It is interesting to note that this 
quantity exhibits a local minimum at q 2 ~ 0.2GeV 2 before dropping to zero at larger 
q 2 . Apparently the frequently employed multipole parametrizations will not accurately 
reproduce this feature of Gs{q 2 )- Kinematical corrections have not been incorporated here 
because they are not relevant for the momentum transfers of the present problem |T[ 
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When SU(3) symmetric wave-functions are employed to compute the spin and flavor parts 
of the matrix elements the strange magnetic moment increases by one order of magnitude 
as compared to the use of the exact eigenfunctions. This is not surprising since the use of 
these wave-functions assumes that strange quarks are as light as the non-strange ones. 
Obviously virtual strange degrees of freedom are easily excitable in that case. 

The realistic amount of flavor symmetry breaking yields a small and negative value 
for the strange magnetic, even smaller than other estimates do [TJ|, [3J|. To some 

extend all these predictions seem to be contradictory to the preliminary results from 
the experiment at MIT-Bates: G S M (0) = 0.46 ± 0.36 ± 0.08 ± 0.18. The first erros is 
statistical, the second due to the background and the third takes account of uncertainties 
in the radiative corrections entering the analysis ||25|| . This measurement apparently 
favors a sizable positive value, however, due to large errors small negative values are 
not ruled out either. On the other hand, a large negative strange magnetic moment, 



as some calculations predict [26, 27], is unlikely. Estimates of the strange axial charge 
G S A (0) are closely connected to the issue of the proton spin puzzle |J. Employing flavor 
symmetry and using data from semi-leptonic hyperon decays results in G S A (0) ~ —0.1. 
When deviations from flavor symmetry are taken into account this number may easily be 
reduced and even a zero value is possible |], [28[ . 



4. Results and discussion 

Before employing the above computed form factors to evaluate R y we will (in fig. ^) 
compare the old energy spectrum for the LAMPF neutrino beam with the new and more 
accurate one ||. As can be observed when comparing the solid and dashed lines the energy 
distribution of the muon-neutrinos have not changed significantly. However, the new 
muon-antineutrino flux is different in shape and magnitude from the old one. Previously 
the latter was estimated by simply assuming it contributed 20% to the total (i.e. v — u+V) 
neutrino beam. The improved spectrum for the antineutrinos is peaking at lower energies 
hence the yield induced by the antineutrinos is significantly decreased. The reason is that 
the extraction of R y requires a threshold of En > 60 MeV for the energy of the knocked 
out nucleon to discriminate against events from elastic neutrino scattering on free protons. 

We have next substituted the data of the improved spectra to predict R y as a function 
of the strange quark axial charge, G S A (0) for different values of the strange quark magnetic 
moment G S M (0). A dipole parametrization, which may not be very realistic cf. fig [I], has 
been employed for the momentum dependence ||. In fig. ^| the results are displayed for 
both the quasi-elastic antineutrino- and neutrino-induced reactions on 12 C. The range 
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assumed for G S M (0) is taken from the estimate of refs. p2| |23|][]- For the non-strange form 
factor also a dipole approximation has been adopted. As the reduction factor is the same 
for the proton and the neutron, the ratio cr(v, v'p) / a(y, v'n) is not effected when using the 
improved spectra. Hence the almost linear dependence on G S A (0), which was previously 
found (c/. fig. 4 in ||]), is recovered, although the absolute antineutrino induced cross 
sections have decreased by a factor of approximately 4.3. Since the interference terms 
contribute with opposite signs to the V- and //-induced reactions, R y exhibits a stronger 
increase with G S A (0) for the V induced reaction than for the one associated with the 
neutrino. 

Fig|| also contains the prediction for R y obtained from the momentum dependent 
form factors predicted by the SU(3) Skyrme model with vector mesons as outlined in 
the preceding section. While, for the case of antineutrino-induced reactions, the ratios 
are slightly different for the three parameters sets discussed above, they are practically 
indistinguishable for neutrino-induced scattering. As a consequence of the quite small 
value for G S M (0) the prediction of this model is expected to be approximately given by 
the dashed lines in fig. |3|. However, this is not quite the case since the model does not 
exactly reproduce the non-strange form factors, cf. table |. 

In fig. |] we finally display R y [y = v + v) as obtained when both the neutrino and 

antineutrino yields are included. As a consequence of using the improved antineutrino 

spectrum, R y {p) is not quite as sensitive to G A (0) as previously || |J estimated. On 

the other hand the sensitivity on G^(0) is increased. Furthermore the negligibly small 

77-flux causes the curves of fig. ^|b and fig. [| to be identical within 1%. From fig. |] we 

deduce 1.1 < Ry{v) < 1.2 as the prediction obtained from the Skyrme model with vector 

mesons. Adopting flavor symmetric wave-functions {i.e. SU(3) D-functions) significantly 

increases this prediction, R y (v) ~ 1.4. Thus an experimental determination of R y {v) will 

in particular serve as a test of the quantization procedure proposed in ref. fl9fl . In this 

context it should be stressed that a precise measurement of R y will especially be suited to 

(dis) approve models for the baryons. The reason is that for this quantity the uncertainties 

inherited from the nuclear model will cancel almost completely. Barbaro et al. |13| have 

recently addressed this issue by comparing the single proton and neutron yields as well 

as their ratio R y in two extreme nuclear models with G A (0) being a parameter. It has 

then been observed that R y is almost independent of the nuclear model but very sensitive 

to G A (0), which, of course, is a prediction of the nucleon model. On the other hand 

the separate yields may significantly vary when G A (0) is kept fixed but different nuclear 

5 More recently a similar analysis which attempts to make contact with perturbative QCD has been 
performed in ref (^4). These authors consider their result G S M (0) = — 0.24± 0.03 as an upper bound (in 
magnitude). However, neither of these calculations incorporates the feature of flavor symmetry breaking. 
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models are employed. 

In our treatment strange degrees of freedom enter solely by virtual excitations of 
strange quark-ant iquark pairs inside the nucleon. Let us also briefly coment on uncer- 
tainties which enter due to other effects. One may e.g. wonder whether or not the coupling 
of the leptonic current to the meson exchange currents between the constituents of the 
nucleus lead to additional contributions. This problem has been addressed by Musolf et 
al. [29 1 for the case of 4 He. Although this nucleus is less complex than 12 C it represents 
a J = and T = nucleus as well. The exchange terms were observed to considerably 
improve on the agreement of the non-strange isoscalar form factor of 4 He with experiment 
(see fig. 2 of ref. ||29|| ). Nevertheless they have turned out not to be relevant for momentum 
transfers less than approximately 500MeV for the strange vector form factor (see fig. 3 
of ref. |^9|). From the energy spectrum displayed in fig. ^| we hence conclude that our 
results will not significantly be modified by the contributions associated with the meson 
exchange currents. From a conceptual point the meson exchange terms are incorporated 
to satisfy current conservation for the many-body system as this requirement is not nec- 
essarily satisfied by the one-body current operator. One may estimate the corresponding 
uncertainties by employing various expressions for the current operator which are related 
by the continuity equation [30|. For different excitations of 12 C Friar and Haxton [30] have 
compared the electric form factors predicted by various forms of the current operator with 
experimental data. They find that the form of the current operator which incorporates 
current conservation constraints only at zero momentum transfer (q = 0) but exhibits 
the correct high q behavior contains at least some meson exchange effects. It is exactly 
this form which enters our calculation. The feature of correct high q behavior is shared 
by neither the form of the current operator which is completely constrained by current 
conservation nor by the one which is not constrained at all. 

We would also like to point out a possible source of uncertainty which is related to 
the measurement of the energy En of the knocked out nucleon, which may be relevant 
for the interpretation of the experiments. Comparing the ratios R y in fig. f| with the 
previous results (fig. 5 in ref. we notice that the curves differ by a constant shift of 
approximately 15% in R y (i)). While for vanishing strangeness contributions to the form 
factors the present calculation (dashed line at G S A (0) = in fig. f|) gives the consistent 
value R y = 1.0, we extract R y ~ 0.85 from fig. 5 in ref. ||. This shift is caused by 
adopting different thresholds for the integrated neutron- yield. In ref. 0] da/dE^ has been 
integrated over the energy of the emitted nucleon starting at En = E p = E n + 2.77 MeV 
> 60 MeV (i.e. as a consequence of the Coulomb-shift the energy threshold of the 
neutrons was taken to be E® = 57.23 MeV), contrary to the results shown in fig. and 
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Table II: The ratios Ry for neutrino-induced nucleon knockout on 12 C (and for G S A (0) = 
G S M (0) = 0) as a function of different thresholds for protons (E®) and neutrons {E®). 



E u p [MeV] 


K [MeV] 


R y {u) 


60.0 


60.0 


1.004 


60.0 


59.0 


0.943 


60.0 


61.0 


1.069 


61.0 


60.0 


0.943 


61.0 


59.0 


0.885 


61.0 


58.0 


0.832 


61.0 


61.0 


1.004 



fig. d, where we have set the thresholds to E N = E® = E® = 60 MeV. This crucial 
dependence of R y on the E p - and .^-thresholds is linked to the steep decrease of the 
neutrino flux (fig. |2|) between 70 and 200 MeV, which also causes the differential cross 
section da/dE^ to strongly decrease with increasing En. We investigated this effect in 
more detail by slightly varying E N for protons and neutrons. The results for neutrino- 
induced reactions are summarized in table [IT] (where we have set G S A = G S M = 0, for 
simplicity). It is apparent that the experimentally determined value for R y crucially 
depends on an accurate measurement of the energy En of the knocked out nucleon. A 
relative shift of only 1 MeV between the thresholds E® and E® will cause an error for 
R y of ~ 6%. We find the same ~ 6% error for non-zero G S A and G S M as well. It should, 
however, be remarked that the absolute value of the E^-threshold has no influence, as 
long as E^ and E® are identical. This independence was already pointed out previously, 
see fig. 3 of Ref. |J. 
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Figure captions 



Figure 1: The momentum dependent strange vector magnetic form factor G S M in the SU(3) 
vector meson Skyrme model (left panel). For completeness the strange axial vector form 
factor G S A is also given (right panel). Different parameter sets have been used, cf. table |. 



Figure 2: Comparison of the previous and updated (anti)neutrino spectra for the LAMPF 
decay-in-flight neutrino source. The solid, dotted and dashed lines refer to the updated 
neutrino, updated antineutrino and previous (anti) neutrino spectra, respectively. The 
antineutrino spectra are always scaled by a factor four. 



Figure 3: Ratio of integrated proton-to-neutron yield for quasi-elastic antineutrino- (up- 
per part) and neutrino-induced (lower part) reactions on 12 C as a function of — G S A (0) for 
different values of G S M (0) within the theoretically estimated regime |22], [23|]. The symbols 
indicate the predictions of the 577(3) Skyrme model with vector mesons, cf. table |. Their 
location on the horizontal axis reflects the associated prediction for G S A (0). 



Figure 4: Same as fig. |3], for the sum of antineutrino- and neutrino-induced yield. 
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